Lower entropy bounds and particle number fluctuations in a Fermi sea 
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In this Letter we demonstrate, in an elementary manner, that given a partition of the single 
particle Hilbert space into orthogonal subspaces, a Fermi sea may be factored into pairs of entangled 
modes, similar to a BCS state. We derive expressions for the entropy and for the particle number 
fluctuations of a subspace of a Fermi sea, at zero and finite temperatures, and relate these by a 
lower bound on the entropy. As an application we investigate analytically and numerically these 
quantities for electrons in the lowest Landau level of a quantum Hall sample. 
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The study of quantum many particle states, when mea- 
surements are only applied to a given subsystem, are at 
the heart of many questions in physics. Examples where 
the entropy of such subsystems is interesting range from 
the quantum mechanical origins of black hole entropy, 
where the existence of an event horizon thermalizes the 
field density matrix inside the black hole [ll l3l j to entan- 
glement structure of spin systems ^. In this work we 
address the relation of entanglement entropy of fermions 
with the fluctuations in the number of fermions. 

First we show that given a subspace of the single parti- 
cle Hilbert space, a Fermi sea may be factorized into pairs 
of entangled modes in and out of the subspace, thereby 
writing the state as a BCS state 15]. 

While upper bounds on entropy, were a subject of 
numerous investigations, especially since Bekenstein's 
bound,4j, lower bounds on entropy are less known. 

We show that given a " Fermi sea" , the entropy of the 
ground state, restricted to a particular subspace A of 
the single particle space, relates to the particle number 
fluctuations in the subspace via the inequality: 



Sa > 41og2AArl > -81og2 < A^l > 



(1) 



Where Sa is the entropy associated with the subspace 
A, and AiV^ are the fluctuations in the particle number 
in A, and <C N\ ^ is the fourth cummulant of particle 
number [l6j. The importance of this result lies in the 
fact that particle fluctuations are, in principle, easier to 
measure, and are fundamentally related to the quantum 
noise in various systems. On the technical side note that 
the right-hand side, has the advantage of being easy to 
calculate analytically in a wider class of problems. 

We start by examining the ground state of non- 
interacting fermions, in arbitrary external potential, 
when measurements are applied to a given part of the 
space. The basic example is a Fermi sea or a Dirac sea 
where we are interested in the relative entropy of a given 
region of space, and in fluctuations in the number of 
particles there, but one may also consider entanglement 
in Fermion traps (Fermi degeneracy of potassium atoms 
C^^K) has been observed by De Marco and Jin 5]). 

The discussion is also relevant for systems which be- 
have like a non-interacting Fermi gas, as in problems of 



transport at the zero temperature limit |l7l|. in ideal met- 
als, where transport may be approximated well within a 
non-interacting theory, due to good screening. In an ideal 
single channel conductor the analogy is done by mapping 
excitations that travel at the Fermi velocity to a time - 
energy coordinate representation in discussion of quan- 
tum pumps la, Li lal ■ The analogy is especially manifest 
in the problem of switching noise 9] . 

The ground state of a noninteracting Fermi gas, con- 
taining N particles is obtained by occupying the allowed 
states (j)i G H {H is the single particle Hilbert space) up 
to energy Ef, i.e.: 






(2) 



Where '0^ sue creation operators which satisfy the usual 
canonical anti-commutation relations (CAR): 






(3) 



The state \gs > or the "Fermi sea", is a typical ground 
state for a large class of Hamiltonians (sometimes it is 
necessary to carry a suitable Bogolubov transformation). 
It can also describe spin chains via the Jordan- Wigner 
transformation. 

Let A be a subspace of the single particle Hilbert space 
H, so that H = A®A^ and let E = span{(f)i; 1 < i < N} 
be the subspace of occupied single particle states of H 
(the Fermi sea). 

Let Pa be the orthogonal projection on A. Consider 
the matrix M(A)y =< Pa(J)j , PA(t>i >, i,j = 1, ..,iv[l3. 
Mij is a hermitian matrix so that it can be diagonalized 
by a unitary U: M — W dia.g{di)U . 

The new orthonormal modes are defined by 



Ai 



Ek ulPA4>k 



where the factor di is the 1-th eigenvalue of M, and serves 
to normalize the Ai with the inner product on H . Sim- 
ilarly we take Bi — ,,"' ,'^ , which are orthonormal 
since M{A^) = I - M{A). Obviously A^ e A and 

B, e A^. 



Since U is unitary we write E = span{W(j)i}- Using 
Uji(j)i = ^fd'iAi + \/\ — diBi, the ground state may be 
written as: 






(4) 



Note that the < d^ < 1 due to the structure of M. In 
the particular case that [Pe,Pa] = then the exercise 
is trivial, as they can be diagonalized simultaneously. In 
addition, whenever there is a symmetry operator that 
commutes with Pa and Pe, then the resultant eigen- 
niodes are invariant under the symmetry. 

The ground state may be written also as a BCS 
state in the following way. Consider the vector \A >= 
nV'^(^i)|0 >, defined by "filling" the modes in A, and 
redefine ipi^i) — V'I(^j)j then: 



\9s>-' 



Ui^ 



^l-d,^HB,)^l{A,))\A> (5) 



Similar to a BCS state, where the entangelment structure 
of the modes pairs A^, Bi is clearly seen. 

If we factor the single particle Hilbert space H into 



N 



H = hR span{Ai, Bi} ff^ (Complement), 



(6) 



the fermion Fock space factors into an appropriate tensor 
product. Wc can write naturally the density matrix as: 



N 



\gs >< gs\ 



Wi >< Vi 



(7) 



where v, = (v^*t(^.) + VT~rfI*^(Sj))|0, >. After 
a partial trace over the B' s, the reduced density matrix 

IS given by: pA = Qs>i^i \ q ^. 

probability of having a particle in mode Ai 
From the density matrix pA we have: 



where d,- is the 



ATVl 



<7V>A=Tr(M), 
=< N"^- < N >^>A= TrAf (1 - M) . 



(8) 



To study further moments, denote P{k) the probability 
of having k fermions in A and consider the generating 
function: 



X(A) 



Z^ 



P{k)e'^'' = dct(l + M(e'^ - 1)). 



(9) 



The cummulants of the number of fermions in A may 
be extracted from the generating function by differenti- 
ating logx and setting A = 0. For example, the fourth 
cummulant is given by: 



«7v4»=a4,iogxU=o = 

Tr(A/(l - M)(l -6M + 6Af2)) 



(10) 




FIG. 1: The functions —a; log a:: — (1 — x) log(l — x) > 
41og2a::(l-a;) > -8 log 2a::(j; - 1)(1 - 6a; + 6a;^) are related to 
the entropy, 2nd and 4th cummulants 



The reduced entropy is obtained from pA'- 

Sa = -Tr(Af log Af+ (1 - A/)log(l - M)). (11) 

This entropy is equivalent to the usual entropy of a Fermi 
gas with occupation number operator M . 

Comparing the expressions l(5|l. (|10|l and Hll() . and using 
the fact that < A/ < 1, we get the inequalities (^ as 
promised (Fig. QJ. 

The difference between the cummulants and the en- 
tropy comes from eigenstates with probability away from 
1/2 but not exactly or 1. The fourth cummulant may 
be used to estimate the contribution to the second mo- 
ments of eigenvalues oiuA away from 1/2, as the function 
—2x{x — 1)(1 — 6a; -I- 6x^) is negative outside the interval 
(1/2 -y/t/A, 1/2+^/3/1) 

The inequalities are valid for finite temperatures too. 
To see this, consider a general density matrix of the form 
p — Z^^e~^'^°>°^. By tracing out the A-^ degrees of 
freedom, the reduced density matrix acquires the form 



PA ~ det(l — nA)e 



log( 



:)• 



(12) 



Where ha = Pajt^Pa is as an operator on A. The 
construction of this distribution from the covariance ma- 
trix < aia, > was discussed in several papers, mainly 
in the context of density matrix renormalization group 
|l3. 14]. where a similar expression was obtained. 

The resulting density matrix on A may be considered 
as thermal by appropriately choosing the " energies" e^ — 
log "^ ' . The density matrix may then be written as a 
"thermal" state p = ^^ig-E ^'V'n^.)V'(A,). Note that 
while the new eigenstates, may have large values for the 
energy with respect to K, the operator ua is well defined, 
and bounded by 1. The entropy is: 

S = -Tr(nA log riA + (1 - "-a) log(l - ha)) 

The generating function x(A) becomes: 

x(A)=det(H-n^(e^^-l)), (13) 

and the expressions for the moments are the same with 
the identification: AI -^ ua- 

One may attempt to find a similar inequality for 
Bosons. For a closed system in thermal equilibrium 

5* = Tr((l -f ube) log(l + ube) - tlbe log ube) (14) 
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FIG. 2: S{R)/AN(R)^ for a disc in the lowest landau level 

(where nsE is a Bose-Einstein occupation operator) and 
{ANf = Ti-nBE{l + nBE)- (15) 

In most cases particle fluctuations will dominate the en- 
tropy, however one can check that whenever hbe < 1 (a 
dilute boson gas) then: 



Sbosous >ilog2)iAN)] 



(16) 



For systems in which the single particle energies have a 
gap A above the zero energy, as the temperature goes to 
zero we have rii — > 0, making the inequality above valid. 

While the cases of spin chains and free fields, were 
studied in numerous works, other important examples 
wait addressing. Here we consider the problem of 2D 
electrons in a quantum Hall sample. 

The filled lowest Landau level states are spanned by 

In [lOJ, the entanglement entropy of a Chern-Simons 
theory describing a quantum hall defined on a disk was 
shown to scale like the radius. Let us now calculate the 
particle fluctuations. We choose A to be a disc of radius 
R. Due to the radial symmetry of the system and of A, 
the factorization I^J is possible mode by mode, and the 
lowest Landau level may be written as: 



\LLL >= l[{^/d;^^PH\k >a) + Vl - dkH\k >aJ)\0 :fl7) 
fe 

Where dk = I \,' — - are given in terms of the in- 
complete gamma function, \k >a= ''h=x{\z\ < R)\k >, 
and I A; >a± is defined similarly |20(. In this case that 
both the entropy and particle fluctuations are asymptot- 
ically linear in R, thereby proportional to the boundary 
area. 

Particle fluctuations in this model are given by 
< AiV^ >= J2i Mii{l - Mil). This sum may be approx- 
imated analytically as follows: write 



EiM',, = j:k,^f^.iMk,i = 



(18) 



^ C de J,"" J,"" dxdye^' 
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Which can be written, after the x,y integrations, as the 
contour integral: 






1 

'27ri 



dz 



(1 



-B^^l-z) 



)(1 



^R^l^l/z) 



\z\ = l 



[i-zy 



FIG. 3: The modes fc = 6n, n = 0, ..7 for _R = 6, 
trated closer to the boundary with increasing k 



are concen- 



Since the integrand is analytic outside z = we break 
it up and do part of the integration on a contour with 

\z\ — !■ oo: 



(A7V(i?))2 = i?2 _ 



27ri 



-"'(e-^ 



-2¥l4| = ld^ 



|2|^00^^ {1-Zy 



{1-zr 



The first integral can be evaluated in the limit i? — > oo, 
by the residue theorem, and simply cancels the R^ term. 
The contribution to the remaining integral comes from 
around the point 9 = 0, we linearize near this point and 
get, as i? — *■ oo: 



{^N{R)r-^lZd9^-^ 



-^R 



(19) 



We see that in this model, the fluctuations are propor- 
tional to the boundary. A numerical check shows that 
the entropy and < AA^^ > are comparable (Fig. EJ. 

In this example it is manifest that for a fixed R, as 
k grows larger (especially |2l|. for k > R^), the modes 
\k >A are localized stronger at the boundary, showing 
that the information is concentrated near the boundary 
(Fig. |3J). Note that in this case we consider from the 
outset an infinite number of fermion modes, in contrast 
with the spin-chain case, where a finite segment yields a 
finite dimensional fermionic Fock space after the Jordan- 
Wigner transformation. 

We remark that since quasi-free fermionic states are re- 
lated to determinantal processes [ll| , equivalent inequal- 
ities may be obtained for the connection of the entropy 
and the fluctuations in such systems. In particular, if one 
is interested in particle fluctuations, then the theorem of 
Lebowitz and Costin |l2l |. ensures a gaussian behavior 
for the scaled particle number fluctuations as the volume 
grows large for a large class of determinantal processes. 

In the case of a ID Fermi gas, both entropy and par- 
ticle fluctuations in a box of size L scale as log kpL |23| 
Particle number fluctuations in ID conductors may be 
observed using ultra- fast transistors, however it will not 
be possible to ignore the Coulomb interaction. We note 
that for a ID conductor, there is a pre-factor due to the 
cross section, which means that the fluctuations may be 
quite large. 



In this Letter we have shown how the ground state 
of fermions can be factored into sets of pairs of modes 
- inside and outside a given subsystem. We have out- 
hned the connection between the entanglement entropy 
and the fluctuations in the particle number in two ways: 
First, the inequality ^ supplies a lower bound on the 
available entanglement entropy, and second, noting that 
in some cases (namely the Id Fermi sea and the lowest 
Landau level) both quantities scale in a similar way. 
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Note that a fourth cummulant may have either sign 
The zero temperature is completely degenerate. Other- 
wise, we need the condition (for strongly degeneracy) 

We restrict ourself to states that obey \Pa4>j > P > 0, 

otherwise they are already of the required form. 

In the symmetric gauge for the vector potential Ai = 



The length unit is V2i where I = J -^ is the 



magnetic length. 

Note that the average number of particles contained in 

yl is X; rffc ~ R^- 

In this problem there is a crossover near k = R^, where 

dfc ~ 1/2. To see this use the property ju^ < 1/2 < 
^<"-"'^\ from which dfl2(i?^) < 1/2 < d^2(i?2-l) (vahd 



r(n) 

for R^ integer) 

when disorder is present we expect the particle fluctua- 
tions to behave like S = log l^kp, where l^ is the elastic 
mean free path. 



